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Abstract 

The Marcus-Lushnikov process is a finite stochastic particle system, in which each particle is entirely 
characterized by its mass. Each pair of particles with masses x and y merges into a single particle at a 
given rate K(x,y). Under certain assumptions, this process converges to the solution to Smoluchowski 
equation, as the number of particles increases to infinity. The Marcus-Lushnikov process gives at each 
time the distribution of masses of the particles present in the system, but does not retain the history 
of formation of the particles. In this paper, we set up a historical analogue of the Marcus-Lushnikov 
process (built according the rules of construction of the usual Markov-Lushnikov process) each time 
giving what we call the historical tree of a particle. The historical tree of a particle present in the 
Marcus-Lushnikov process at a given time t encodes information about the times and masses of the 
coagulation events that have formed that particle. We prove a law of large numbers for the empirical 
distribution of such historical trees. The limit is a natural measure on trees which is constructed from a 
solution to Smoluchowski coagulation equation. 

1 Presentation of the problem 
1.1 Introduction 

Let E = (0, oo). Let K : E x E — > (0, oo) be a symmetric continuous function. Let 5* be 
the set of finite integer- valued measures on (0, oo). S contains elements of the form 



x = 



for n G N where yx, . . . , y n > are distinct and for i G {1, . . . , n}, m,i G N. The Marcus- 
Lushnikov process with coagulation kernel K is the continuous time Markov chain on S 
with non-zero transition rates given by 



q(x, x') 



m i m j K(y i ,y j ) if i < j 

hniimi - l)K(y h yi) if % = j 



fOr X' =X + 5y i+yj -5 yi -5y.. 

Let us give a way of constructing a Marcus-Lushnikov process (X t )t>o- Let [N] = 
{!,..., N} and let yi, . . . , y^ > be the masses (not necessarily distinct) associated to 
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each particle in [N] . Set 

N 

i=i 

For each % < j take an independent random variable such that is exponential with 
parameter K(yi,yj), and define 

T = minTy. 

i<j 

Set X t = X Q for t < T and 

X T = X t - (5y t + 5y 3 - Sy i+yj ) 

if T = Tij, then begin the construction afresh from X T . Each pair of clusters i < j 
coagulates at rate K(y h yj). 

Now fix hq a measure on E and take a sequence (xq) n>q E S such that 

= N-'x" - (1) 

weakly on E as A" — > oo. Let (X^) t > be Marcus-Lushnikov with kernel starting from 
x^. For each N, we can write 

m(N) 
i=l 

with m{N) < N and , . . . , > non necessarily distinct. Set 

^ = N-'X t N . 

Without loss of generality, for the rest of the paper, we will take m(N) = N. Indeed, 
after one step in the process above, the number of masses in the system will be less than 
iV and so we will be exactly in the case where m(N) < N. 

Our aim in this paper is to set up a historical analogue of the process and to prove 
that it converges to a limit measure that can be constructed from the strong solution to 
a generalised form of Smoluchowski's equation [l],[2](to be made precise below). 

Before defining precisely this new process let us explain why it is interesting to know 
about the history of formation of a cluster. 

The Marcus-Lushnikov process [5] describes the stochastic Markov evolution of a 
finite system of coalescing particles. It gives at each time the distribution in masses 
of the particles present in the system but does not retain any other information that 
the particles might contain. In other words, we lose in part the information contained 
in the particles that is their history. Why is it interesting to know about the history 
? For instance, consider a system of N particles with associated masses yi, . . . , yjv > 0. 
Assume that these particles can only be of three types say either A, B or C. Allow them to 
coagulate according to the rules of coagulation of the Marcus-Lushnikov process. Then, 
the usual Marcus-Lushnikov will give us at each time the masses of the particles present 
in the system but will not be able to tell us for each particle present at this time how 
many particles of type A, B or C this particle contains along with the order of formation. 
Our historical measure will give us at each time the particles formed with their respective 
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masses, the time when they have formed but also the history of the formation from its 
beginning that is the time at which each intermediary particles have formed along with 
what they contain. We could think of an other application in industry in the process 
of making a certain chemical product. We can assume that in order to make a certain 
powder we need to put N ingredients in a specific order and at specific times. Then our 
historical measure will allow us to follow the formation of the powder and to detect if 
ingredients were put in the wrong order at the wrong time. 

We are now going to review the work of [1] and [2] about the convergence of fj,f to 
the strong solution to Smoluchowski's equation as we will use this tool to prove our main 
result (stated in 1.3). 



1.2 Related work 

Take ip : E — > (0, oo) to be a continuous sublinear function. We suppose that ip > 1. 
Assume that the coagulation kernel can be written as follows: 

K(x,y) = K(x,y)<p(x)<p(y) (2) 

with K bounded on E x E. For // a non-negative Borel measure on E such that, 

/ K(x,y)n(dx)n(dy) < oo 

J ExE 

we define L(/i) as follows: 

(/, L(jm)) = \( {f(x + y)- f(x) - f(y))K(x, y)iM(dx)n(dy) 

A J ExE 

for all / bounded and measurable. We consider the following measure-valued form of the 
Smoluchowski coagulation equation, 



Ik = l*o + 



I L(fi s )ds. (3) 
Jo 

We admit as a strong local solution any map: 

t^nt: [0,T)^M + 
where T e (0, oo] and A4 + is the set of non-negative Borel measures on E, such that 
1. for all B C E compact, 

t^Ht(B) : [0,T) -> (0,oo) 

is measurable. 

2. 

/ / (p 2 (x)fi s (dx) < oo 
Jo Je 

for all t < T. 
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3. for all bounded measurable functions /, for all t < T, 
/ f{x)nt{dx) = / f(x)fio(dx) 

JE J E 

+ (f(x + y) - f(x) - f(y))K(x,y)fx s (dx)fx s (dy)ds. 

JO JExE 

In the case T = oo we call a strong local solution a strong solution. Assume that 

(y? 2 ,/i ) = / (p 2 (x)/j, (dx) < oo. (4) 

JE 

Then [2] tells us that there exists a unique maximal strong solution to (3) denoted (fJ,t)t<T 
for some T > 0. Moreover assume that 

/ f(x)<p(x)no (dx) -> / f(x)(p(x)n (dx) (5) 

JE JE 

as A" — > oo for all / bounded and continuous on i?. Then, for all t < T, for all / 
continuous and bounded on E 

I f(x)<p(x)fi?(dx) -> / f(x)(p(x)nt(dx) 

JE JE 

as A" — > oo in probability, that is 

AT 

as A" — > oo weakly in probability with weight function </?. Indeed, for e > 0, for all / 
continuous and bounded on E, for all t < T, we can find v > 0, such that for all A/" > 0, 



P ( sup 



f(x)<p(x)fi 8 (dx) - / f(x)ip(x)n s (dx) 

E JE 



> e J < exp(-i/7V). (6) 



1.3 Our main result 

Assumptions 

The following assumptions hold for the whole paper. Let /j,q be a measure on (0, oo) 
and take a sequence (Ho) N>0 (as defined in subsection 1.1) such that (1) is satisfied that 
is 



1 N 

= n S S v" ^° weakl y as A^ oo. 



fc=i 

.A 7 



For convenience we will write ^ = y& for all /c e [AT] = {1, . . . , A^}. Let E = (0, oo). 
Let K : £ x E — > (0, oo) be a symmetric continuous function. Let : E — > (0, oo) be a 
continuous sublinear function. We suppose that tp > 1. Assume that (2), (4) and (5) are 
satisfied and that K is bounded on sets of the form [0, M] 2 with M > 0. Let (^t)t<T be 
the strong solution to (3) with T > 0. 
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Basic notations for Trees 



For a finite set J C N, we write T(J) for the smallest set containing J, such that 
{i,j} G T(J) whenever i,j G T(J). We refer to elements of T(J) as trees. They are finite 
binary trees with leaves labeled by elements of J. Each i G T( J) has a set of leaf labels 
A(i) C J determined by 

A(i) = {*} 

for i G J and 

A({i,j}) = A(i)UA(j) 

for all i,j G T(J). For each % G T(J), |A(i)| will denote the number of elements of the set 
A(i). Let n : T(J) — > N be the counting function, defined as follows: 

n(i) = 1 

for i G J and 

for all i,j G T(J). Finally, for J C N finite, we define 

T*(J) = TV = {i G T(J) : |A(i)| = n(i)}. 
This is the set of rooted binary trees with leaves labeled by distinct elements of J. 

The coagulation process on Trees 

According to the rules of construction of the usual Marcus-Lushnikov process, we set 
up an analogue of this process, each time giving the trees present in the system. We 
start with the set of labeled particles [N] = {1, . . . , N} with masses y 1 , . . . , > 0. For 
i = {^1^2} £ T*[iV], the mass yi of % is given recursively by 

V% = Vh +Vi 2 - 

Each tree particle i G T*[iV] has a tree type r{i) G T, where T = T({1}), determined 
by t(i) — 1 for % G [N], and r({i,j}) = {r(?),r(j)} for i,j G T*[iV]. We construct the 
historical analogue of the Marcus-Lushnikov process on T*[N] as follows. Set I[0] = [Nj. 
For % G [N], set = 0. There exist trees . . . ,i[N — 1] G T*[JV], coagulation times 
< ■ ■ ■ < Si\N-i] < oo (they are exactly the jump times in the usual Marcus-Lushnikov 
process), and subsets . . . , I[N — 1] C T*[iV], such that for each n G {1, . . . , N — 1}, 
(A(i) : % G I[n\) is a partition of [TV], and I[n] satisfies the recursive relation 

I[n] = I[n - l]\{ii[n],i 2 [n]} U {i[n]} 

where = ^[n]}. 

Let us give an example with N = 5. The drawing below represents one of the config- 
urations we can obtain from the coagulation process on trees described above. 
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Si[l] Si[ 2 ] Si[ 3 ] Si[4] 

Here, 7[1] = {{1, 2}, 3, 4, 5}, 7[2] = {{{1, 2}, 3}, 4, 5}, 7[3] = {{{1, 2}, 3}, {4, 5}} 
and 7[4] = {{{{1,2},3},{4,5}}}. 

The historical measures 

Fix t < T. Our principal object of interest is a process of empirical particle measures 
jif on the space of historical trees A(0, t) which we shall now define. The space A(0,t) 
is given by 

A(0,t)=\jA T (0,t) 
where Ai(0, t) = (0, 00) and for r = {ti, t 2 } G T, 

A T (o,t) = {(s, {6,6}) : s e (0,t),6 g A n (o,s),6 e A T2 (o, s )}. 

Let us illustrate these historical spaces through an example. Take r = {{1}, 1} G T. 
The tree below represents a £ = (s 2 ,{6,6}) £ A-(0,i) with 6 G ^4{i}(0, s 2 ), 6 £ 
Ai(0,s 2 ). 




t time 



Here 6 = (s 1: {y u y 2 }) , 6 = 2/3 and f = ^s 2 , {(«i, {yi, y 2 }), ys}] ■ 

We equip A(0,t) with its Borel cr-algebra (we explain in Appendix 5.1 how to equip 
A(0,t) with a topology). We define on A(0,t) the mass function m : A(0,t) — > (0, 00). 
For £ G Ai(0,*) = (0,oo), we set 

™(0 = 6 

Recursively for r = {ti,t 2 } G T, for £ = (s, {6)6}) £ ^r(0,t), we set 

m (f) = m(6) +m(6)- 
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The empirical historical measure fif is given by 



iei(t) 

where I(t) C T*[N] is the set of trees present in the system at time t. I(t) is given by 

I(t) = I[n\ for S i[n] < t < S i[n+1] . 

For % G [N],$ = yi and for % = {h,i 2 } G I(t), with $ = s, we set Q = (s, {C,C 2 })- As 
we trace back the past of a particle we obtain a tree. Observe that this empirical measure 
jif and our usual Marcus Lushnikov process jif (defined in subsection 1.1) are related 
through the following equality, 

N ~N -1 

We are interested in taking the limit of this empirical measure as iV — > oo. We define 
the limit measure on A(0,t) as follows. For £ e Ai(0,t), we set 

P*(dQ = exp ^- J K(y,y')n r (dy')dr^ Ho(d£) 

where y = m(£) and {^ r )r<T is the strong deterministic solution to (3). Recursively for 
r = { Tl , r 2 } G T, £ = (s, {6,6}) e ^r(0, *) with s < t < T, we define 

P*(dQ = e(r)i^(m(£i),m(£ 2 ))/i s (rf6)/i s (^2)exp ^- J K(y,y')fi r (dy')dr^ ds 
where y = m(£), e(r) = 1 if t\ ^ r 2 and e(r) = | if ri = r 2 . 
Our main result 

Our aim in this paper is to prove the following result. 
Theorem 1.1. For all t < T, 

^ - in (7) 

weakly on A(0, t) in probability. 

This theorem is proved is Section 4. Before giving in subsection 1.5, an outline of 
how we are going to prove Theorem 1.1, we need first to introduce some more material. 
The next subsection 1.4 is dedicated to introduce a coupled family of processes built 
on the same probability space JF, P) (to be specified below), which we shall see are 
Marcus-Lushnikov and which will be really convenient to use for most of our intermediary 
proofs. 

1.4 A coupled family of Marcus-Lushnikov processes 

We start with the set of particles [N] with associated masses yi, . . . , yjv > (that is /Iq). 
For J C [N] we set T* J = T*J\J. 
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A probability space and some random variables 

Let Q = (0, oo) T + [JV] . A probability measure P on Q is denned for u = : % G T*.[N]) 
by 

F(du) = || exp(— uo,j)duJi. 

ieT* [AT] 

Let (Ui : i G T+[7V]) be a family of random variables on Q denned as follows. For 
u = (uji-.ie T* [N]) G tt, set U^uj) = Wi . Then, under P, (Ui : % G T*.[N]) is a family of 
independent exponential random variables with parameter 1. Set JF = cr(Ui : « G TFj_[jV]). 
Let us now define on Q a family of random variables from which we will construct our 
coupled family of Marcus-Lushnikov processes. For i G [N] set 

3 = 0. 

For {i,j} G T*[iV], define recursively 

N 

S {iJ} = max(Si, Sj) + ^(y- y-) U{iJ} ' 
The coupled family of processes 

We are going to build a coupled family of processes ^(X/) <>0 : J C [jV]^ as follows. 

Fix J C [iV] and set nj = |J|. Set J^j = cr(£/; : i G T\J) and fi J = (0, oo) T + J . Fix 
uj — (ut : i G TPj. J) G £V. We start off with the set of particles J% — J with respective 
masses (yi : % G J). We set 

and we consider 

w° = min 

with i 7^ j 

Almost surely, we have v° = S{i j }(uj) for some unique io,jo G J% with io 7^ Jo- We then 
obtain a new set of particles = J^\{io, jo} U {{«0j jo}}- 
We set 

X t J (u) = X{ 

for t < v and 

Now, starting from J^, we consider 

= min 

with i ji j 

Almost surely, we have v 1 = Sj^j^u) for some unique G with ^ 7^ j^. We 

obtain a new set of particles J% = J^\{ii,ji} U {{«i, Ji}}- We set 

x? (u) = x^) 
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for v < t < v 1 and 



ieJ% 



We start again as above with J% and so on. The process stops when there is only one 
particle left in the system. 

Therefore, for each u G VL J , there exist trees . . . ,i[nj — l] G T*J, coagulation times 
Si[i] < ■ ■ ■ < Si[ nj -i] < oo (they are some of the Si(u)'s), and subsets J^, . . . , J^/^ 1 C 
TV, such that for each n G {1, . . . , nj — 1}, (X(i) : i G JJJ-) is a partition of J, and J% 
satisfies the recursive relation 

JN = JZT 1 \{ii[n],i2[n]}U{i[n]} 
where i[n] = i 2 [n]}. We set 

{Si[ n ] if i G J^T^Jn 
oo ifi = i[nj-l] 
SM if^U^- 1 ^ 

T/ can be interpreted as the time of death of the tree particle i in J, if this particle was 
alive in the system at some time. If J = [N], we set t\ n ^ = T { . The empirical historical 
measure on trees fif can be rewritten as follows 



^"jy E 1 {Si<t<T i }S^i. 



i€T*[N] 

The Marcus Lushnikov property for this family 

Theorem 1.2. Let J C [JV]. TTie process (X/) zs Marcus-Lushnikov with kernel ^ 
starting from 

In particular, ^xj^j and (X/ v ) <>0 Ziaue same distribution. 
Proof of Theorem 1.2 : Starting from 

= E ^ ' 

it is clear that the first jump has the correct distribution for Marcus-Lushnikov. Let us 
now look at the fc + jump. We condition on 

4k = ( *o = £ *i = E ■ ■ > x i = E 6 * 

where for I G {1, . . . , k}, J l N is the set of particles present on [si, sj+i) and si < • • • < 
are the jump times. What are the transition rates to go from step k to step k + 1? 
Conditional on we consider the set of particles and look at 

T' = min <% j} = min s {hj} + V {hj} 

with i ^ j with i ^ j 



where for each i,j G J% with i ^ j, V^jy is exponential with parameter K ^ y ^ and, 

' if i,j e J . 

Si if % was formed at some Sj G {si, . . . , s^} and j G J. 

'{m} = ^ s i ^ 3 was f° rme( i a ^ some Sj G {si, . . . , and i G J. 

max(sj, Sj) if i was formed at some Sj G {si, . . . , and 
j was formed at some Sj G {si, . . . , Sfc}\{sj}. 

Thus, for each i,j G J% with i 7^ j, V^jy has started at time s^jy and has been running 
for a duration of Sk — Nevertheless, by the memoryless property for exponential 

random variables, for h > 0, 

P + Vfij} > S k + h\s {iJ y + V {iJ y > s k ) = P (V {iJ y > h) . 

Thus, it is equivalent to add — S{jj} to Sj^} and consider that V^jy starts from s^. 
Hence, conditional on Ajt, we find exactly the T we considered in the construction of the 
the Marcus-Lushnikov process in 1.1. Therefore, the transition rates are the same. Thus 
(X/) t>Q is Marcus-Lushnikov. 

□ 

For all t > 0, we set = N~ l X^ . Note that by construction, for all J C 

[iV], (/if r ' J )t>o is measurable with respect to Tj = cr(Ui : 2 G T* J). We shall see that the 

coupled family of Marcus-Lushnikov processes ^(W^) : J ^ will be useful in 

most of the intermediary proofs leading to our main result. Let us now give an outline 
of the intermediary results we need in order to prove Theorem 1.1. 

1.5 Outline proof of main result 

Fix t G [0, T). In order to prove Theorem 1.1, we need to prove that for all / G Cb (A(0, t)) 

(f,tf) = [ fm?(do - (z./i*) = / mwo (8) 

JA(0,t) JA(0,t) 

as iV — > 00 in probability. To prove (8), we shall see that it is sufficient to prove that for 
all r G T, for all / G C b {A T (0,t)) 

= / fmfrn - </,&> = [ fmm (9) 

JA T (0,t) JA T (0,t) 

as N — > 00 in probability. Then we will be able to conclude, using a tightness argument 
(that is explained in subsection 4.1). To show (9), it is sufficient to prove that 

E ((/, ft »-</, ft) (10) 

and 

E((/,/^) 2 )^(/,/i t ) 2 (11) 
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as N — > oo. Then, 



E((/,/if-/i t > 2 )^0 



and a fortiori 




in order to obtain (9). Finally in Section 4 we will use the results out obtained in the 
previous sections to prove Theorem 1.1. 

2 Intermediary computations 

For the whole section we start with the set of particles [N] with associated masses 
yi,...,yjv > (that is and we consider the probability space P) that we 

have defined in subsection 1.4. We fix < t < T. The aim of this section is to calculate 
E «/,/2f» for / G C b (A T (0,t)) with r G T. 

2.1 A finite sum of conditional expectations 

Let us fix r G T with n leaves that is n(r) = n and take / G Cb(A T (0, t)). We are going to 
express E ((/, p,^)) as the expectation of a finite sum of conditional expectations. Recall 
that for % G T*[iV], r{%) denotes the type of the labeled tree i. We can write, 



as / is supported on A T (0,t). For J C [JV], and fceN define 

T* J = {i G TV : n(i) = k}. 

For any i G T*[iV] with type r, there are 2 q( - T ^ permutations possible of the particles 
composing % which will keep the tree % invariant where q(r) is the number of symme- 
tries in the tree r. This is given recursively by q(l) = and for r = {ti,t 2 } G T, 
q(r) = q(n) + q(r 2 ) + l{ Tl =r 2 }- 

Example : Take % = {1,{2,3}}. This tree has for type r = {1,{1}} and q(r) = 1. 
The permutations leaving the tree i alike are the identity permutation and the one send- 
ing 1 to itself and exchanging 2 and 3. Hence 2 <?(r ) permutations leave this tree invariant. 

Now fix i G T*[n] with type t(iq) = t. Define 




iST* [JV] 





ieT*[JV] 

with r(i) - 



= r 



S n ,N = \ injections a : [n] — > [N] 
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For a G S n> N, set <r(l, . . . , n) = (<x(l), . . . , cr(n)) . Let 

[N]\ n] = {(^, . . . ,i n ) G [N] n :h,...,i n distinct}. 

Observe that for each (ii, . . . , i n ) G W]* n ] there is a unique a G S^at such that c(l, . . . , n) = 
(ii,...,i n )- Hence, 

[NY [n] = {a(l,...,n):aeS n>N }. 

For a G S^jv define o"(i ) G T*[iV] to be the tree obtained from i by replacing each 
particle i G [n] in i by <r(i). 

Example : li N — 5, i = {1, {2, 3}} and a G S 3:5 is such that er(l) = 3, cr(2) = 5, cr(3) = 1, 
then 

<7(i ) = {3, {5,1}}. 

Define 

Jo : [iV]f n] - T*[7V] by I (<r(l, ...,n)) = a(i ). 

Thus, we can write 

2— <?(t) , . 

{S CT (i o) <t<T CT(io) } 



N 

{h,...,i n )e[N]* n] 



Hence, setting J(cr(i )) = [iV]\{A(cr(io))}, we can write 



2-<?( T ) 



We shall compute 



E(/(^ 0) )l{5 CT(i0) < t <T CT(i0)} |^ 0)) ) (12) 

for all a - G S^jv- In order to compute this quantity, we shall find it useful to work with 
labeled trees. 

2.2 Working with labeled trees 

We are going to introduce spaces similar to A(0,t), but for labeled trees. 
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Vector tree of masses 



Each particle j G [N] has a given mass yj. Take a tree i G T*[AT]. Assume that this 
tree is formed from the particles %\ < ■ ■ ■ < i n . We write y = (y^, . . . ,Vi n ) for its asso- 
ciated vector of masses. Its vector tree of masses y is the tree of masses obtained from 
i by replacing each particle in this tree by its mass. For instance i = {1, {2,3}} and 
V = (yi, 2/2, ys) give y = {yi, {y 2 , y 3 }}. 

Labeled historical spaces of trees 

Define Af (0, t) for k G {1, ... , N}, by 

Af(0,t) = {k}x{y k } 
and recursively for % = {11,12} G T£[iV] define 

A^t) = { (s,{6,6}) :«G(0,tUiG 4(M,6 € <(0, S )} 

where y,?/ 1 and y 2 are the respective vectors of masses of and ii- For i G T*[iV] 
(without any associated masses) define 

A(0,t)= |J A?(0,t) 

ye(0,oo) n 

and set 

00 

A(P,t)= |J |J A(0,t). 

N=l ieT*[AT] 

Observe that when we integrate over A^(0,t), we only integrate over the coagulation 
times as the masses of the particles are fixed whereas integrating over Ai(0,t) means 
integrating over the coagulation times along with the masses. 

Reduction of the problem 

Each particle j G [N] has a given mass yj. Take r G T with ti(t) — n for some n G N. 
Fix / G Cfe(A T (0, £)). In order to compute E ((/, /i^)), without loss of generality, by 
subsection 2.1, we need to compute 

e(/ (e*) i {Si < t<Ti} i^ j ) iv- 1 

for i G T*[n] with type r(i) = r and associated vector of masses y = (yi, . . . ,y n ) and 
J = [N]\[n]. The map 

9i : A^t) ^ A T (0,t) 
on forgetting labels is 2 9<T ) to I. Define 

fi = f ° 9i- 
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We have ft G C b (A^O, t)). Let Q e A T (0, t). Then, 

g7 1 (Q) = {CteM0,t):g i (Ct) = Q}. 

The set g^iQ) contains 2 9<T ) elements. Hence, we can write 
So, 

E(/K)l{«,<t<31})=2- 9(T) E( £ /(<*(Ct))l{*<*<7»)- 

If we fix (® G g~ l {il)-, then all the other elements of g^iQ) can be obtained from £ t ° 
by permuting masses between the particles that form a symmetry in the tree i. Hence 
<?j -1 (£t) can be written as a set depending only on (£. 

Example : Take % = {1,2}. Then r(i) = t = {1}. Let Q = (s,{y,y'}) e A T (0,t). 
Then, 

9;\%) = {(*, {{1} x {y}, {2} x {y'}}), (s, {{1} x {y>}, {2} x {</}})}. 

It is clear in the set g^iQ) that we can obtain one particle from the other by exchanging 
masses between particles 1 and 2 which are symmetric in this tree. 

Therefore, 

E(/ ($ !{*<*<*!>) = (Ct°) 1{S,<*<21})- 

Hence, 

e(e(/ ($ I^tjI^ 7 )) = e(e(/< (C t °) I^^toI^))- (13) 

For convenience, according to the context and the spaces we consider, Q will stand for 
either an element of A T (0,t) or Ai(0,t) and more particularly since the masses of the 
particles are fixed as an element of Af(Q : t) where y is the vector of masses of the tree i. 
Hence it is enough to compute 

e (ft (ei) L {Si<t<Ti} ■ (14) 

We are going to compute (14) first for some particular i and ft. Then, we will use these 
intermediary results to solve the general case i G T*[ra] and ft G C& (Ai(0,t)) and by the 
relation (13) we will obtain an expression for f £ Cb (A T (0,t)). 

2.3 Case / = 1 and r = 1 

We take r = 1 G T and / = 1 G Cb (A-(0, t)) . The corresponding set of labeled trees with 
type 1 is [N]. Without loss of generality take % — 1. The corresponding ft G Cb (A(0, t)) 
defined in subsection 2.2 is j\ = 1. Also, since % G [N], we have by definition Si = 0. 
Hence, in this case, we want to compute P(7\ > t | JFjJ where J\ — {2, . . . , N}. 



14 



Theorem 2.1. Let J ± = {2, . . . , N} and T Jx = a (Ui : i G T* Ji) . Then, 

P(T\ > i | JjJ = exp ^- jf J^K( yi ,y)fi^(dy)dr^ a.s 

where 7\ «s as defined previously the time at which particle 1 dies in (/i^) r<t and y 1 is 
the mass of particle 1. 

Proof of Theorem 2.1 : Since /if r ' Jl is measurable with respect to JF Ji; we can write 
^f> Jl = F/^E/i : i G T* Ji) where Jf : (0, oo) T + Jl -> (0, oo) is a measurable function. 
Fix u Jl = (w/ 1 : i G T* Ji) G (0, oo) T + Jl and set fi?' Jl (v Jl ) = Jf (w/ 1 : % G T^Ji) (this 
notation will be kept for the whole paper). 

There exist ri(c<j Jl ) < ■ ■ • < r n {uj Jl ) such that ji^ Jl {uj Jl ) is constant on the interval 
[rk(uJ Jl ), rk+i(u Jl )) for k — 0, . . . , n with the convention r (u Jl ) = and r n+ i(u Jl ) = t. 

Moreover, we can write 

Ti = Gf ((C/ {M} : i G TVi), {Ui : % G T;Ji)) 
where : (0, oo) T * Jl x (0, oo) T + Jl — > R is a measurable function. Hence, 

r 1 (w Jl ) = Gf((?/ {lii} :ierj 1 ),w Jl ) 

is a random variable. To obtain F(T 1 > t | J-jJ, by Fubini, it is enough to compute 
P(Ti(cj Ji ) > t). 

Now observe that 

{Ti(cu Jl ) >t} = {Ti(co Jl ) > t} n {Ti(co Jl ) > r„(cu Jl )} n . . . D {Ti(u; Jl ) > n(uj J ^} 
because 

{Ti(c Jl ) > t} C {Ti(c Jl ) > r„(c Jl )} C . . . C {Ti(c Jl ) > n(a; Jl )}. 

Hence, 

n 

P (Ti(c Jl ) > i) = JJP (Ti(c Jl ) > r fe+ i(^)|Ti(^) > r fc (^ Jl )) . (15) 

Let us start by computing the quantity P(Ti(cj Ji ) > ri(c<j Jl )) . Then, we will calculate 
P(Ti(^) >r fc+ i(^)|Ti(^) >r k (cu^)) for k G {1, . . . , n}. 

For {i,j} G T*[iV] set 

AT 
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Observe that {Ti(cj Ji ) > n(uj Jl )} = {S {1J} (cu Jl ) > ri(u; Jl ) : j = 2, . . . , N} and that for 
j G {2, . . . , N}, S {hj} (u J i) = V {1J} . Hence, 

P (T!^ Jl ) > =P > n(cu Jl ) : ./ 2 V) 

j=2 ^ ' 

= exp \-\^ N r ^ ) J 

by independence of the (V{i,jy '■ j G Ji). Now, for all r G [0,ri(w Jl )), 

3=2 

So, 

^ ri(wJl) ^(yi,y)^ Jl (^)(dy)dr- = £ 

Hence , 

P(T!(o; Jl ) >n(^)) =exp } j^iyi, y)^' Jl (u Jl )(dy)dr j . 

Now, let us compute P (T\(a; Jl ) > r fe+1 (cu Jl ) | Ti(u; Jl ) > r k (uj Jl )) for fc G {1, . . . , n}. 
For k G {l,...,n}, let J k (oo Jl ) be the set of particles present in [^ yh (uo Jl )) for r G 
[rfe(cJ Jl ), rfe + i(cJ Jl )) and assume that j k is the particle formed at time r k (uj Jl ). We can 
write j k = {j'fc(l), Jfc(2)} with j k (l), j k (2) G J fe _i(^ Jl ). Observe that S jk (u Jl ) = r k (u Jl ). 
Then, 

P (T^ 1 ) > rfc+x^- 71 )!^^- 71 ) > r k {u- h )) 

= P (5 {1J} (o; Jl ) > r k+1 (u' h ) for all j G J^)^^) > r fe (cu Jl )) . 

Observe that J k (u Jl ) = {j k } U J k -i(u Jl )\{j k (l), j k {2)}. By the memoryless property for 
exponential random variables, 

P (TUcc^) > r fe+1 (u; Jl ) | T^) > r k {u- h )) 
= P (S {1J} (u' h ) > r k+1 (u' h ) for all j G J fe (cu Jl )\{j fe }, 

5 {1Jfc} (u; Jl ) > r k+1 (uj- h ) | %, } (co Jl ) > r k (u' h ) for all j G J k (^)\{ Jk }) . 

Now we can write J k {uJ Jl )\{j k } = I k {uj Jl ) U K k (uJ Jl ) where K k (u Jl ) C {ji, . . . ,jk-i} 
and I k {oo Jl ) C {2, ...,iV}. For j G I k {uJ Jl ), we have S^ij^a/ 1 ) = and for 

j G ^(cu" 71 ), there exist r fej G {ri (a; Jl ),..., r k _i{uj Jl )} such that S , j (c<j Ji ) = r kj . So, 
S'{i >i }(cj Jl ) = Sj(u Jl ) + V{i tj y = r kj + V{i tj y. Hence, for each j G J fe (^ Jl )\{j fe }, only 
depends on V {1J y. Also, Sj fc (^ Jl ) = r k (u Jl ) so S , { i iifc} (cj Jl ) = r fc (w Jl ) + V {1Jk y. Thus, 
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using the independence of (V{±j},j E Jk(u Jl )\{jk}) and the memoryless property for 
exponential random variables, we obtain that : 



P(T 1 (W /1 ) > r k+1 {u- h ) | T^) > r k {u' h )) 

= ^(V {ljk} >T k+1 {u- h )-T k {u- h )) 

n p ( v i^} > r ^ ji ) i vat} > r ^ ji )) 

II P fau} > r k + i(uJ Jl ) - r kj | V {1J} > r k (uj J i) - r k .) 

jeK k (wh) 

II V (V { i, j} > r k+1 (u J i) - r k (u J i)) 
II P(V{i J }>r fe+1 (cc; Jl )-r fe (cu Jl )). 



Hence, 

P(7\(w Jl ) > r fe+1 (cu Jl ) | T\(cu Jl ) > r fc (cu Jl )) 

= n p(^{u}>^ + i(^ ji )-^ ji )) 



jeJ fe (ur 7 i) 



Now, for all r G [r fc (cu Jl ), r fc+ i(cu Jl )), 



J] exp (-^M (ri+1 (^) - ))) . 



r +liWl) I K(y 1 ,y)tf'*(u>*)(dy)dr = £ (r fc+1 (^) - r fe (^)) . 



So, 

T fc+1 (w J l) 

'r fc (« J i) 

Hence, 

P(T 1 (a; Jl ) > r fc+1 (co Jl ) | T^o/ 1 ) > r fc (cu Jl )) 

r.r fc+ i(o; J i) /• \ 

/ K( yi ,y)tf>*(u*)(dy)dr . 

r fc (w J l) •/£ / 

and using the equality (15) we obtain : 

P(T 1 (o; Jl ) > t) = exp (- j ^K( yi ,y)^ J ^)(dy)dr^ . 
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= exp 



Finally, by Fubini, we get : 

P(7\ > t | T Jx ) = exp (- J* J E K( yi ,y)^(dy)dr^ . 

as required. 

□ 

2.4 Case / and r general 

Recall that each particle j G [N] has a given mass i/j > 0. Take r G T with n(r) = n, 
for some n G N and fix 2 G T*[n] with type r(i) = r and associated vector of masses 
y = (yi, . . . , y n ). Fix / G Q, (A r (0, £)) and consider /j G 0,(^(0, £)) constructed from / 
as in subsection 2.2. Let J = [N]\[n] and set Tj = cr(Uj : j G T+J). Our aim in this 
subsection is to compute 

where Tj, the time of death of the labeled tree i, has been defined in the subsection 1.4. 
The notation "E^jv.j" means that the expectation is taken conditional on starting with 
the set of particles [N] with respective masses yi,...,yjv > 0. This condition can be 
rewritten as follows : y = (y u ...,y n ) and = ± J^jeJ S vr 

2.4.1 Useful notations 

Fix < t < T. 

The space A(f) 

For k G [N], for £ G (0,t), define 

A(fl = (0, *)*{£}. 

Recursively for i G T*[ra] with associated vector of masses y, for £ = (s, {£i,£ 2 }) £ 
AV(0,*), define 

A(0=A(6)UA(£ 2 )U(MU{£}). 
The projection map A(£) — > (0, i) 

For fc G [iV], for £ G Af (0,t), define the projection map n = IT 5 : A(f) -> (0,t) by 
for s G (0,t), 

n< ((*,{£})) = a. 

Recursively for i G T*[n] with associated vector of masses y, for £ = (s, {£1,^2}) £ 
Af(0,t), define the projection map IT = lie : A(£) — > (0,t) by 

f %H if«eA(£i) 
n(«) = n^u) = { n 6 ( M ) if u eA(&) 

[ r if « = (r,0 e [s,t)U{£} 
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The projection map from A(£) onto the space of trees of masses 

For k G [N], for £ G Af(0,t), set y t (g) = y k and define 

y = y(0 : A(f) -> (0, oo) by y^)(u) = y u = y k . 

Recursively for i G T*[n] with associated vector of masses y, for £ = (s, {656}) G 
^(0,t), set 

w(0 = {i/.(6),y.(6)} 

and define 

1/ = y(0 : A(£) -> T(0, 00) 
where the space T(0, 00) is given by 

T(0,oo) = [J T T (0, 00 

reT 

where Ti(0, 00) = (0, 00) and for r = {t±, r 2 } G T, 

T T (0,oo) = {y = {y\y 2 } : y x G T n (0,oo),y 2 G T T2 (0,oo)}, 

by 

f y(Si)(u) if«eA(eO 
2/(0(«) = y(6)(«) if«e A(6) 

I 2/ u = M&),y s (6)} if «= (r,0 G [s,t)U{e> 

The notation 

For £ G A k (0,t) with A; G [iV], set 

K/: = 1. 

Recursively for i G T*[iV], for £ = (s, {6,6}) G Ai(0,t), define 

2.4.2 An expression for the conditional expectation 

Set T\ n \ = &(Uj : j G TPj_[n]). The aim is to compute 
\^' J (/< te*) 1 {fifi<*<r i }|^\7 V .F [n] ) 

= fi (Ct) 1 { Si <t<TM} E y,l£"' ( l {t<Ti}\rj V Fin]) 

because fi (Q) l{ S . <t<T Wy is measurable with respect to T[ n \ and {t < T { } C {t < T/™'} 
(t/™' is defined in 1.4). The lemma below will be proved later on. 

Lemma 2.2. With the notations above, on the event {Si < t < Tj^}, 

P n,j (t < Ti\Tj V F [n] ) = exp ( - / / K(y r , y')^ J (r (dy')dr ) a.s. 
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Hence, integrating over A?(0,t) we obtain 

ifi ($ 1{*<*<^}I^» = / fi (0 Py < * < if ] , Q e df) 

^AV(0,t) V 7 

6XP (~ Ja(o X y '^n,?(r)( d y') dr ^) • 

Therefore, we need to compute P w ^ < t < T\ n \ Q G d^j for some £ G (0, t). 
Lemma 2.3. 

where v = uf is the Lebesgue measure defined in the Appendix 5.2 and for £ G Af(0,i) 
and s G (0, £), 

r/en^(s) 

mi/i IT^ : A(£) — > (0,£) as defined in 2.4-1- 

Proof of Lemma 2.3: Let us fix £ G Af(0,t). To £ we can uniquely associate 
{(si,ji), . . . , (s n -i, jn-i)} where si < • • • < s n -i are the coagulation times and ji, . . . ,j n -i 
are the labeled trees (subtrees of i) formed at s±, . . . , s n _i respectively. Denote by J k (C,) 
the set of trees (subtrees of i) present on [s fc , s^+i) for G {0, ... ,n — 1} with the con- 
vention s = and s n = t. For each k G {1, ... ,n — 1}, we write j k = (j fc (l), j fc (2)} 
with j fc (l),j fc (2) G J^). Observe that J fe (£) = U J fc ' 1 (0\{j fe (l),J fc (2)}. For 
convenience, for each k G {1, . . . , n — 1} we will write 

= K (yjk(i),yjkw)- 

We start with 



jv,[n] _iv-c 

- N 2^°vr 



j'=i 

Let (Yk)k=i n-i be the jump chain associated to the Markov-Lushnikov process iyh'^ J 

Let Ji, . . . , J n _i be the jump times. Observe that after the n — 1 th jump there is only 
one particle left in the system, and so the process is in an absorbing state. Thus, the 
configuration £ can be represented through the jump times and the jump states by, for 
k — 0, . . . , n — 1, 

For k G {0, ... , n-2}, the rate of going from F fc = ± Ej e jfc(£) ^ to ^fc+i = jj Ej G jfc+i(£) K 

is q k ,k+i = ~ J w L and the rate of leaving the state Y k = ± Eje S yj is <?fc = | E j,j e J fc (5) ■ 

with j ^ I 
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Hence, 



P. 



, ,(5 4 <t<3f , ,€{€de) 

Ji G . . . , J n _i G Y\ = — ^ 5 yj , . . . , Y n -\ = — ^ <5 % J 



P„ 



ds 



= qo,l ■ ■ • 9n-2,n-l ex P ■ ■ ■ ex P (-5n-2(s„-i ~ S„- 2 )) d«l 

Replacing and gfc^+i by their respective values we obtain that 

/ \ 

i K (yi>yi) 



n-l- 



P w (Si<t<TW,$edt) 



TT — 

11 /V 



fc=l 



cxp 



V 



j'.ieJ 1 ^) 

with j I 



N 



Si 



I 



( 



... x exp 



1 

2 



\ 



iV' 



n-l 



exp 



with j ^ i 

1 m) 



N 



Sn-l — Sn-2) 







N 



ds uf(dO 



as required. 



□ 



Now let us prove Lemma 2.2. Before proving it, let us introduce some basic no- 
tations. For oo = (ujj : j G Tj_[iV]) G Q, we will write a;' 71 ' = (a;j : j G T+[n]) and 



or 



= jGT*J). 



Proof of Lemma 2.2: Fix£ G A^(0,t). We can find o^"' = : j G T* [ra]) representing 
the configuration £. To £ we can uniquely associate {(si,ji), • • • , (s n -i, jn-i)} where 
si < • • • < s n -i are the coagulation times and ji, . . . ,j n -i are the labeled trees (subtrees 
of i) formed at Si,...,s n _i respectively. Denote by J fc (£) the set of particles from £ 
present on [s k , s k +i) for k G {0, . . . , n — 1} with the convention s = and s n = t. For each 
fc G {1, ... ,n - 1}, we write j k = {j k (l),j k (2)} with j k (l),j k (2) G J^" 1 ^)- Observe that 
J k (0 = {h} U J k jfc(2)}. Also, for convenience, for each A; G {1, ... ,n — 1} 

we will write 

K ik = K (yj k (i),yj k m)- 

Now, fix u£ = (uf : i G T* J). There exist n(uj$) < • • • < 7v_i(u;^) such that ^' J (u^) 
constant on [r k (u)Q) , r k+ i(u>Q)) for k — 0, . . . , n' — 1 with the convention ^(cUq) = and 
r n i{uJo) = t. For each G {0, ... ,n' — 1} denote by Jfc(o'o) the set of particles from 
/if ,J (cj^) present on [r k (u)$),r k+1 (u$)). Define 

Q(o; J ,0 = G : cj w = u%\uj j = 4)- 

To obtain our result, conditional on {Si <t< xf 1 '} it is enough to compute : 

F y,^ J H w G fi (^o , : * < TM)) = P w> ^ (* < ^(4, 4" 1 )) (16) 
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for short where T^uj^ujq) is understood as a random variable on Q(uq,£). Since we 
are conditioning on {Si < t < T^}, on each interval [sfc, s fc+i)> the particles present at 
time Sk in £ that is J k {!;), will not interfere between each other until time s k +i- Hence, 
conditional on {Si <t< T^}, the quantity (16) means that on each interval [sfc,Sfc+i), 
each particle in J fc (£) will not interfere with the configuration ^/i^' j (uJq )) r>s until s^+i- 
This is equivalent to consider on [s fc , s*i+i), for each i G J h (t;), starting from 

Jfc(s fc )U{i } _ A 
ieJ fc (s fc )u{j()} 

that i wm survive until s k+ i where Jk(s k ) is the set of particles present in (this 
set is one of the J p {oOq) for p G {0, . . . , n' — 1}). If jl^ Sfc ) u ^ ^ denotes the time of death 
of the particle i when starting from the set of particles Jk(s k ) U {i } (this death time 
is defined in 1.4) we want T l J o k{sk)u{io} (u) > s k+1 for all u G It is exactly the 

result we have obtained in Theorem 2.1. For short we will write, j'^ k ^ Sk ' )U ^ ^ — . Thus, 
setting 

A 1 (co J ,0 = {cuE fi(o; J ,0 : V* G J° '(£) , 7f H > «i} 

A n _!(a; J ,0 = G n(atf,0 : V* G J"~ 2 (0, 7?-» > 
A„(cu J ,0 = {cu G ^K J ,0 : 7? n _» > *}, 

we can write, conditional on {Si <t< T/™'}, 

(* < T *KV N )) = (Ai(wo J ,0 n . . . n A n (u J ,0) ■ 

The process we consider is a continuous-time Markov process whose state space is the 
space of measures on the set (0, oo). Hence, fixing the coagulation times s± < • • • < s n _i 
we can consider the process independently on each interval [s k , Sfc+i) starting respectively 
from (y k , ^ where y° = (y 1 ,...,y n ) and for k G {0, . . . , n - 1}, y k = (yj,j G 
J k (0). Thus, 

p w ^ (* < T i(" J o^o ] )) = F y,^Hl) (M<*U)) x P ,V^) ( A =»K J .0) x • • ■ 

■■■ x Vv&W) K^o.O)- 

First, let us compute 

p ^'ho (AK J ,0) = p ({^ g n( Wo J ,o : v* G J (0,^» > si}) . 

Take i G J°(£). For j G Jfc(^o) for k e I 1 ; • • • X) either = V{ io j} if j is a 

simple particle, or there exist p G {0, . . . , k} such that Sjfafi) = t p {oJq) if j is a composed 
particle. In this case >% ,j} = t p {uJq) + ^Wi}- Hence, when i G {1, . . . ,n}, T^(o>q , u;^) 
only depends on (\/ {loj} : j G J k {uj^, k = 0, . . . , ra'-l). Hence, (t/(cu j , 4" 1 ) : J e AO) 
are independent random variables. Thus, 

^W^^OH n p »^)H wGn ( w »'^ : W >ai H- 
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For j e J°(0, by Theorem 2.1, 

PrfW^ G : T ^ U) > Sl}) = 6XP I - /„ I^K{y 3 ,y')^> J H)(dy')dr j . 

Hence, 



'0 J E 



P iV 



Now let us compute P fc jv.j, j., (^(cug , 0) f° r ^ ^ {2, . . . , n — 1}. We are starting at 

time Sjt_i with the set of particles J fc_1 (£)- By the memoryless property for exponentially 
distributed random variables, it is equivalent to start off at time with the set J fe_1 (£) 
and to go to s k — s k -i- Hence, the computation of P t , n,j. iAA k (uji, £)) is the same than 

for F y ,^' J (u,J) (^i( w o>0) replacing J°(£) by J fe_1 (0, ji by j fc , and si by s k -s k -i. Thus, 



P 
Thus, 



v»tf> J M)(M«>o, 0) = exp f - f* / K( yj M> J (uj J )(dy')dr 



fc=i \ Jsk - 1 JE jeJ k -HO 



Hence, 



= ex P E-/ / E K (vvy')^r ,J H){dy')dr 
= exp (~ /,( ? ) X x ^ r ' V '^n(l)( d y') dr 

P iv, J (t<T^ J V^ [n] )=exp(- / / K{y r ,y')y% J , r) {dy')dr j a.s 
V Ja (® Je St J 



as required. 



□ 



Hence, 

= L ( o/ i(0 ^ exp {~ L®f E KM) ^ r)W)dr ) exp (rl^ ds ) vm) 

= [ m)K,h N ({^ j ) r<v t)vm) 
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where 



ds 



hN ((^' J ) r<t , e) = exp (- J j E K(y r , y')^ J {r) (dy')dr^ exp (- jT ^ 
and where for y G (0, oo) n , the Lebesgue measure v\ is defined in the Appendix 5.2. 
2.4.3 Convergence for teh conditional expectation 

Fix y = (jji, . . . ,y n ) G (0, oo) ra . This is the vector of masses associated to the particles 
{1, . . . , n}. Let y be the tree vector of masses associated to i G T* [n\. Define 

PtMv) = Ptfi(y) = [ fi (0 K € exp (- f [ K(y r , y')^(r)(dy')dr) ^(df) 



where 



Set 



h (Mr<t > = eX P ( ~ / / K (Vr, V')^U 

\ Ja(£) Je 



(r) 



(dy')dr 



Theorem 2.4. 



p t N fi(y) -Pt/i(?/) m probability as N -> oo. 



Proof of Theorem 2.4 ■ For a given £ G (0, t), we can write 

/ \ 

_i ^ K (yj,yi) 



n-l 



^((^• j ) r<t ,e)=ri ex p 



fc=0 



n-l 



2 ^ iV 

i i,ieJ fc (0 

\ with j ^ £ 



■(s fc+ i - S fc ) 



/ 



{J exp 



fc=0 



£ K( yj ,y')^> J (dy')dr 
Sk JE jeJHO 



where si < • • • < s n -i are the coagulation times associated to £ and ji, . . . , jn-\ are its 
labeled subtrees (subtrees of i) formed at these respective times. J fc (£) represents the set 
of tree particles (subtrees of i) from £ present on [s fc , s^+i) for A; = 0, . . . , n — 1 with the 
convention s = and s n = t. As iV — > oo, for each fc G {1, . . . , n — 1}, 



n-l 



/ 



O exp 



fc=0 



\ 



2 ^ AT 

\ with j ^ / 



(Sjfc+l - Sfc) 



1. 



/ 



Now we want to show that 
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n-1 



JJexp | - 

fc=0 



Sk+1 



s * JE jeJHt) 



K(y 3 ,y')^ J (dy')dr 



n-l 

n 

fc=0 



cxp - 



'fc+1 



K(yj,y')fi r (dy')dr 



as iV — > oo in probability. It is sufficient to show that for each k G {1, . . . ,n — 1} and 

(r s k+i r \ / /•sfc+i /• 

- / / i/)tf> J (dy')dr - exp - / / K{ yj ,y')ii r {dy')di 
J s k JE J \ J s k JE 

as iV — > oo in probability. For fc G {1, . . . ,n — 1} and j G J fe (£) consider the random 
variable 



fSk + l f s k+l 

Xtj>= / K( yj ,y')^ J (dy')dr- \ K( yj ,y')^ r (dy')dr. 

Js k 



s k 



On [s k , Sfc+i), 



Hence , 



, at , n,j i 1 x 



/ K( yj ,y')^(dy')dr 
rsk+i 

J Sk 



K(yj,y')n r (dy')dr - (s k+1 - s k ) ^ K ^vVi) 

ieJ k (0 



Now, we know that 

K (yj>y) = K(yj,y)<p(yj)<p(y)- 

For y G (0, oo) define f yj : (0, oo) — > (0, oo) such that 

fvjiv) = K(yj,y)v(yj)- 

By hypothesis, this function is continuous and bounded. Hence, we can write 

K (yj,yi) 



vk,j 



s fc+l 

Sk 



Sk + 1 
Sk 



(fy 3 ¥,Vr) dr - / (f yj ¥,Vr)dr - (s fc+ i - s k ) 



ieJ fc (0 



AT 



Hence , 



TV 



< 



(Sfc+l - s fc ) Y 



K (yj,yi) 



ieJ k (Z) 



N 



+ (s fc+ i - s fc ) sup I (/ % </?, /i^) - (/ %¥ >, /i r ) I 
»*e[s fc ,s fe+ i) 
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The first term converges to 0. By (6), the second converges to in probability as iV — > oo. 
Hence — > in probability as iV — > oo. Hence, for all £ G Af(0,t), 



^((^• J ) r<4 ,e)-fc((^) r<t ^) 



in probability as AT — > oo. 

Now, on ^4|(0, i) _£Q is bounded, say by some C > 0. Thus, 



AV(0,t) 



A? 



\P t N fi(y)-P t fi(y)\ <C\\M 
Moreover, for all £ G AV(0, t), 

Hence by the Bounded Convergence Theorem, 



< 1. 



E 



(7 ^((^ J ) r<t ^)-fc((Mr) r < t ,e) 

\-MV(0,t) V - / V - / 



as AT 
as N 
as AT 



oo. Thus, 



oo and a fortiori 



E(|P t ^(y)-^(y)|)-o 



p?Mv) - PtMv) 



oo in probability as required. 



□ 



3 Some convergence results 

The aim of this section is to prove (9, that is, for all r G T, for all / G C b (A T (0,t)), 

as A" — > oo in probability. 

3.1 Convergence of the expectation 

We recall that we are working with the set of initial particles [N] = {1, . . . , A^} with 
associated masses yi,...,yjv > 0. Let r G T with n leaves, that is n(r) = n. Take 
% G with type r. Associate the map / : [N]^ — > T*[AT] as defined in 2.1. For 

(ii, . . . , in) G [A^]* n j with associated masses y^, . . . , yi n > we will write y^,. ..,«„) for the 
vector of masses and y^,...,^) for the vector tree of masses associated to J(ii, . . . , i n ). Let 
us fix / G C b (A T (0, t)). We aim to prove the following result. 
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Proposition 3.1. 

E({f,P?))^(f,iH) 

as N — > oo. 

Proof of Proposition 3.1 : For M > 0, define : (0, oo) — > (0, oo) to be 

f 1 if x < M 

m M (x) = I -x + (M + l) if M < x < M ■ 1 
[ if x > M + 1 

We can write, 

(/,/if> = (f^ M (m),^) + (/(l - ^M(m)),/if ) 
where m : A(0, £) — > (0, oo) is the mass function defined in subsection 1.3. Hence, 

E ((/, /if)) = E {(f* M (m), /if » + E - * M ( m) ), /if )) . (17) 
We use the following lemma that we shall prove below. 
Lemma 3.2. 

E((/* lf (m) ) jif»^(/* lf M 1 /i t ) 

as N — > oo. 

Thus, taking the limsup and liminf over N in the expression (17) and applying Lemma 
3.2, we obtain, 



Urn sup E «/,/if» = (f* M {m),ik) + limsup E ((/(l - * M (m)),tf)) (18) 

N N 

liminfE((/,/if>) = (/* M (m),^)+liminfE((/(l-* M (m)),/if)). (19) 



AT TV 

and 

lii 

Since /if = /if o m^ 1 and </? > 1, we can write 



|(/(l-^ M (m)),/if)| = |(/(l-M/ M ),/if)| 

<||/||oo|Mi-*m),^>| 

Hence 

E (|(/(1 - * M (m)),ji»)\) < II/IUE (|M1 - * M ),//f>|) • (20) 

Now, I (</?(! ~~ ^m), A*f ) | is bounded for all iV. Indeed, since (p is sublinear, for all t, we 
have 

Hence, 

(V9(l-^ M ),/if)<(^^)- 

Now fix e > 0. Since (</?, /xf) — > (</?, A*o) as iV — > oo, we can find iVo > such that for all 
iV>iV , 
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Let M = maxAr<jv (if, fi^) and set Co = max(M , (<p,Ho) + e). For all N we have, 

{<P,V?)<C 

and a fortiori, 

(^(l-*Af),^><C 

so it is bounded for all N. Moreover, it converges to \(<p(l — \I/m), /•*<)! in probability as 
N — > oo. Thus, by the Bounded convergence Theorem, 

E (| Ml - - E(\(<p(l - 

= \(<p(l-V M ),nt)\ 

as N — > oo. Taking the limsup over iV in the expression (20) we obtain 

Urn sup E (|(/(1 - * M (m)),/if)|) < H/IU - *m),M| . 
iv 

As M — > oo, y?(l — ^m) \ and is positive. Hence, by the Monotone Convergence 
Theorem, 

\(<p(l-V M ),Ht)\ ^OasM^oo. 

So 

limlimsupE (|(/(1 - * M (m)),jj? )|) = limliminf E (|(/(1 - * M (m)),ji?)\)=0 
Moreover, 

|(/(1 - *M(m)),/i t )| < H/IU |M1 - V u ),m)\ 
that is going to as M — > oo. Thus, 

lim (fV M (m),ji t ) = (f,fit). 

Hence, taking the limit as M — > oo in the relations (18) and (19) we obtain 
limsupE((/,/if>) =liminfE((/,/if)) = (/,&). 

AT iV 

Thus, 

as required. 

□ 

Proof of Lemma 3.2 : Set / M = f^/ M (m). We can write, 

»«/".#» E e(/«(^<« «)i„ n , „,<«,„ „„) 

(ii,...,«n)e[JV]f B] v 7 
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where I : [N]* n ^ — > T*[iV] is the map defined in subsection 2.1. For (ii, . . . ,i n ) € [-W]* n ]> 
let be the corresponding element for f M in Cb(A I ^ li ___ tin )(0 : t)) as defined in 

subsection 2.2. Then, 



E«/ M ,/^» 
1 

" iV 



E E 

(h,...,i n )e[N]* n] 



fM ( >J(ii,...,i„)\ -i \ 

/J(ii,..,i„) ^ J i {S /(ll ,..., ln )<t<T /(lli ...,, n) } I 



(ii,...,i„)e[JV]f B] L^'ln .J' '*' 



exp - 



A(0 



^.^^^-■^(dlOdr) exp (~ f^ hA ^ n) > 



using the expression of the conditional expectation obtained at the end of the subsection 
2.4.2. 



Step 1: We are going to give an alternative expression of E ((/ M ,/2^)) using a new 
measure. Define on (0, oo) n the measure, 



n,N 



- E 



(n,...,i„)e[iV]* 



For (ix, ...,?„)€ [iV]f n] define 
G(yn,---,y in ) 



E 



%,-,>n) ^ 7 



AT 



ds 

(dy')dr ) ^^(dO 



In this expression, the first three terms in the integral (inside the expectation) only 
depends on (y^, . . . , y in ), whereas the last term depend on (y 1 , . . . , y N ). Nevertheless we 
only need to know which particles are present in the first three terms to work out which 
particles are going to be in the last term. Moreover, if there exist G [N] with %\ ^ i[ 
but = y^ then we have 

G{yi' v yi 2 ...,y in ) = G{y h ,y i2 . . . ,y in ). 

Thus we can write, 



1 



G( yil ,...,y in ) 

(i 1 ,...,i n )e[N]* n] 
/ G(x u . . . ,x n )nl ,N (dxi . . . dx n ). 

J(0,oo) n 
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Step 2: We are going to show that we can write 

E «/ M tf)) = I G(x 1 , . . . , x n )^(d Xl ) . . . tf(dx n ) + o(l). 

J(0,oo)« iV 

Let us consider 

/ n y 
where Cq ,jv is a finite sum of terms of the form 



1 N 

N,m,...,n k _ L \^ r 



ii,...,ifc=l 
with ii, . . . , ij. 
distinct 



times times n k times 



for fc > 1 with ni + n2 + . . . + nfc = n modulo some permutations of the masses. Take 
fc < n — 1 and consider, 



N 



1 v / ^i"(ii,...,n,...,ifc,...,ifc)^ 

ii,...,ifc=l 



»l,...,Zfc 

with ii, — , i fc 
distinct 

1 {57( n ,..., n ,..., lfc ,..., lfc) <t<T /(ili ... iili ... ilfei ... ilfc) }^^ r 

where J(ii, . . . , i ± , . . . , i k , . . . , i k ) is obtained from ...,i n ) with (h, . . . , i n ) G [N]* n] 
by substituting in the tree . . . , i n ), the particles {i p : 1 < p < rii} by the particle i±, 
the particles {i p : ri\ + 1 < p < ri\ + n 2 } by the particle i 2 - ■ ■ , and finally the particles 
{ip : n — n k + 1 < p < n} by the particle i k . It is a tree where particles are not all distinct. 
We are going to prove that this quantity is o(-^). By a similar argument to before we can 
write 

N 



h,...,i k =l 
with ii, . . . , i k 
distinct 



where g : (0, oo) fe 
Now define 



is the expectation above. 



W) 



N 



i=i 



and on (0, oo) fc the measure, 



N 



k,N 



^ ^ ^{Vil i---iVi k 



h,...,i k =l 
with ii, . . . ,i k 
distinct 



Then, 



ti N < MO 
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We can identify ^' ni, - ,nfc and /j,q ,N . Hence, 

f g{x l ,...x k )$ N (dx 1 ...dx k ) 

J(0,oo) k 



tt N N r, 



- l^k / 9(xi, ■ ■ ■ x k )i$ (dxi) ...fiQ (dx k ) 

iV J(0,oo) k 



For i e T*[n\, let A itM (0,t) = {£ e Ai(0,t) : m(£) < M + 1}. On A iM (0,t), the map 
A(0, i) — > (0,oo) : £ — > is bounded by C = (sup J/)y ; g [ 0)M ]2 2/'))" 1 ( which is 
attained as K is bounded on [0, M] 2 compact). Outside A,m(0, t), ^uijn) is the zero- 
function and so is /m- Now, for all (ii, . . . ,i n ) G [iV] n non necessarily distinct, 



Now, 



g(y h , . . .,y ik ) = E (P t N 'f% u ... tin) (j/(u, ...,;„))) ■ 

\ P t N fl(i u ...,i n ) {y(h,-,in)) I 

-14 ,><„,/" »®««»(-f,zf*) 

cxp (-LX K <^^><*^) , fcw<*) 

■ /A /(il,...,i„)V U ' t J 

= C||/l|oo / 5^ l{ a i<...< a „_i}dsi...ds n _i 



(0,oo) Possible permutations 
for si, . . . , s„_i 

< CH/IU"- 1 . 



Hence, 



J(0,oo) k 



7V " -/(0,oo) 



Cll/ll f 1-1 /" 

N i(0,oo)fe 

_ c\\f\U-\ N k 

- Nn _ k VP, ) ■ 

as ip > 1. It is clear that this quantity is going to as N — > 00 since < n — 1 and 
(<P,I1q) — > ((p,fJ>o) < 00 as iV — > 00. Hence, 



E((/M,/if» = / G(x 1 ,...,a; re ) / i^(^ 1 ).../i Af (^ ri ) + o(l). 

7 (0,00)" iV 

as required. 
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Step 3: Let i e T* [n] with type r and without associated masses. For each (xi, . . . , x n ) E 
(0, oo) n , we can associate to iq the vector of masses (xi, . . . , x n ). Hence we can write 

G(x 1 ,...,x n )=E{P»fg(x 1 ,...,x n )) 

By Theorem 2.4, for all x±, . . . , x n e (0, oo), 

P t N f%(x u ...,x n )^P t f%(x u ...,x n ) 

in probability as A" — > oo. Also, P^f^(xi, . . . ,x n ) is bounded by C||/|| 00 t n_1 . Hence, 
by the Bounded Convergence Theorem, 

E «/m, /if)) -> / ^/if ((&i, • • • , a;„)) M^i) • • • M dx n) 

J(0,oo) n 

as A" — > oo. 

Now, we need to prove that 

/ p t/? (0&i, • • • , z„)) A*o(<fei) • • • A*o(<fcn) = (/ M , fit)- 

«/(0,oo) n 

But, 

/ fa, ...,x n ) Hoidxi) . . . fi (dx n ) 

= [ I ) K^^)^- Xn \d^{d Xl )...^{dx n ) 

exp ( - / / K(y r ,y')n r (dy')dr) 

= / (0 = / / M (0 = (/ M , /i*> 

^A i() (0,t) -/A T (0,t) 

using Appendix 5.3. Hence, 

E((/ M ,/2f))-(/ M ,/i 4 ) 

as A^ — > oo as required. 

□ 



3.2 Convergence of the expectation of the square 

We are still working with [N] = {1, . . . , A^} with associated masses y±, . . . , > 0. Take 
r G T with n leaves that is n(r) = n for some n > 0. Fix / e C&(^4 T (0, i)). The aim in 
this section is to prove the following result. 

Proposition 3.3. 

M</>/^> 2 H</^> 2 

as A^ — > oo. 

In this aim, we are going to proceed similarly to sections 2 and 3. First we are going 
to express E ((/, flf) 2 ) as the expected value of a finite sum of conditional expectations, 
then we will compute these conditional expectations and give some properties about their 
convergence. Finally we will prove Proposition 3.3. 
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3.2.1 An expression for the expectation 

Take i G T*[n] with type r and associate the function / : [N]^ — * T*[iV] defined in 2.1. 
The aim is to compute 

E((/,^) 2 )- 

We can write, 

2~9( T ) 



E«/,/ir> 2 )=E^^ r e / (^) w t <T„ (1 



AT2 



0> 



+ E (^ W ) f W ) 1 {^(0<*<r„(0} 1 {Sa'(0^< T «r'(i)}) • 

<x,<x'€S'„,jV 
with Jm(<r) n Im(cr') = 

where Im(a) = : j G Observe that in the formula above we omitted terms of 

the form 

e^/ ($r w ) / (^r w ) i{5 CT(i) <t<T CT(i) }i{5 CT , (i) <t<T CT , (i) }^ 

for a - , u' G S^jv with Im(a) r\Im(cr') ^ 0. Indeed, in this case, by the way we constructed 
the {Sj : j G T+[iV]}, the events {S a ^ < t < T a ^} and {S a /^ < t < T a i^} are disjoint 
(because by construction the trees have all distinct leaves) and so the quantity above is 
equal to 0. Let us look at the first term in the expression of E ((/, fif) 2 ). We have 



2 -2(?(t) 



AT2 

°"Gi> nj jv 

From subsection 3.1, as TV — > oo, 

E((/ 2 ,^))-(/ 2 ,^)<oo. 

Hence as N — > oo, 

2 -2g(r) 



AT2 

For cr, o"' G S^jv with Im(a) fl Im(a') = 0, set 

J((7(i),(/(<)) = [JV]\ {A WO), A ((/(<))}. 

Hence, 

E ((f,tf) 2 ) = e(^ E e (/ (^ W ) / 

with ira(cr) n Im(a') = 



1 {S CT(l) <t<T CTW }l{5 CT , w <t<T CT , w }|^ J(,TW ' ,T ' W) ^ + 0( — ) 
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Thus, we need to compute the following conditional expectation 

E(/ (C W ) / l { 5 CT(l) < t <T CT(l)} l { 5 CT , w < t <T CT , w} |^ (CTW -' W) 

for a, a' G S HtN with Im(a) H Im(a') = 0. 

3.2.2 An expression for the conditional expectation 

In order to compute E ((/, ft,^) 2 ) , without loss of generality, by subsection 3.2.1, we need 
to compute 

for ^ G T*[n],2 2 G T*[2n]\[n] both with type r. The quantities and y 2 represent 
the respective vector of masses (as defined in 2.2) for %\ and i 2 , J = {2n + 1, . . . ,N} 
and Tj = a (Uj : j G T* J) . Define = / o g h and / ia = / o ^ 2 where ^ and g h 
are the map on forgetting labels defined in subsection 2.2. Then fa G (A il (Q : t)) and 
4 eC 6 (A i2 (0,t)). 

Step 1 : A symmetry argument 

By a similar symmetry argument that the one we used in subsection 2.2, we can write 

E yi 1 y2,^ J (/ (& n ) / (£ 2 ) 1 {5 il <«T il }l{5i 2 <t<T i2 }|^ 

Thus it is enough to compute 

Step 2 : Simplifying the computation 

Set T\ n \ — a (Uj : j G T^[n]) and .F[2n]\[n] — ° (Uj ■ j G T^_[2n]\[n]) and look at 

E y\y 2 ,^-' (/ii fe^ /<a (^) 1 {Si 1 <t<T il }l{5 i2 <t<T i2 }|^V V .F[ n ] V .F[2n]\[n]) • (21) 

Observe that (£ t n ) 1^ <t<r M} is measurable with respect to T\ n \ and that 

/«2 (Ct 2 ) 1{5. < t <T [2nlMn1 } * s measurable with respect to T\2n\\\n\- Moreover, {t < T^} C {t < 

T^ ] } and {* < T i2 } C {t < T% nNn] }. Thus, the quantity (21) is equal to 

fil (Ct 1 ) fii (Ct 2 ) 1 {S^-CKT^} 1 {S l2 <t<Tll n]X[n] } 

E y\y*,^' J (h^T^itKT^lfj V T {n] V .F[2n]\[n]) ■ 
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Hence we can write, 

E y\y 2 ,^* J {fii fe' 1 ) /*2 fe 12 ) l{S il <t<T il }l{S ia <t<T ia }|J : j) 

= / , 2 fh (6) /* (6) 

p w i,^ (X < * < ^ir 1 , 5 fa < t < rif ]Mn] , e e g . 

Lemma 3.4. 1. Conditional on the event {S h < t < T^ ] } n {S i2 < t < t|" 1M " ] } ; 




where 

K s {n\m = -\ E ^(^)- 

Pyi,* < * < T i"U 2 < t < Tt^iZl 1 e ^i, G dfa) 

= -^Mr ex p ^- y o -^v - rfs J ex p J o ~^r ds ) < ^< ( d &- 

In the lemma above, the quantity (1) is similar to the one in Lemma 2.2 except that 

1 2 

here, it means that for £i e (0,t),^2 G (0,t) and c<j j = (a;/ : i G T+J), the config- 
urations £i and £2 do not interfere with each other and are not killed by (/ i ^' J ) r>0 ( ujJ )- 
Hence, mimicking the proof of Lemma 2.2, we obtain the formula of the lemma. The 
second part of the proof comes from applying the Lemma 2.3 along with an argument of 
independence. 
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Hence, 

E y\y 2 ,^' J {fii &) fa &) 1 <*<^ x > 1 {^ 2 <*<^ 2 } I ^) 



f 



(0,t)xAV ( ,t) W .A) iV 

exp (~ X c ) Se K ^ Vri y '^^i (r) ( dy '^ dr ) 

exp (- / / K(y r ,y')^ {r) (dy')dr) v( (d^ (d&. 

V J A(^) J E / 

3.2.3 Convergence for the conditional expectation 

We keep the same notations that in the previous section and we define, 

PffMfiaiV*) = E y 1 ,y 2 ,Ho' J (#) h (ff) hs^T^hs^T^ j) N 2 ^ 

and 

Pthiy^fiAy 2 ) 

= 1 1 2 ^Jn(6)exp (- / K{y r ,y')^ J {r) {dy')dr) uf (d^) 



K bfi2 (6) exp ^- jT ^ y X (y r , y')^n ( J 2 (r) u t ( d &)- 

Proposition 3.5. i. i^/uG/ 1 )/^ 2 ) -> PtfiAy^fhiy 2 ) as N -> oo m probability. 
2. PtfiMMy 2 ) = Ptfniy'WUy 2 )- 

We are not proving this proposition since the proof is similar to the proof of Theorem 

2.4. 

3.2.4 Proof of Proposition 3.3 

We are going to sketch the proof of this proposition as it is very similar to the proof of 
Proposition 3.1. 

Step 1: For M > 0, define ^ M ■ (0, oo) -> (0, oo) to be 

( 1 if x < M 

V M (x) = < -x + (M + l) if M < x < M ■ 1 
[ if a; > M + 1 

We can write, 

= (/^M(m),/if) + (/(l - * M (m)),/^>- 
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Hence, 

(/, = (f^ M (m),^) 2 + (/(I " *mM), tf? + 2</*M(m), $ > (/(l - <Mm)), /2? ) 
and so 

E «/,/if ) 2 ) = E {(FH M (m),ii? ) 2 ) + E - ^(m)),/if ) 2 ) 
+ 2E ((f* M (m),ji» - * M (m)),/if )) . 

Step 2: We prove that E ((/^ M (m), /if ) 2 ) -> (f^ M (m), fi t ) 2 as iV -> oo. Let n 6T*[ti] 
and « 2 £ T*[2n]\[n] both with type r, and associate the maps I\ : [-N]* n ] — > T*[7V] and 
I 2 : [^]f 2n]XW -»• T* n [N] defined in subsection 2.1. For M > 0, let / M = f* M (m) 
where m is the mass function. Let f?f and / 4 f be the correspondent representant for /ir- 
respectively in C b (A il (0,t)) and C b (A i2 (0,t)). Define on (0, oo) 2 ™, the measure 

1 N 

2n,N _ 1 \ ^ r 

~~ _/y2n {yn,-,yi„,yh ,-,%„) ■ 



»l,-,>i>,il,-Jn=l 

with ii, . . . , j„, ji,. . . , j„ 
distinct 



We can write 

E«fVf> 2 ) 



= E ( / P/Vif ((&1, • • • , Zn)) ((l n +l, • • • , &2n)) iff'" (dx U dx 2n ) J 

\^(0,oo) 2 " / 



(0,oo) 
+ °( — ) 

Hence, by mimicking Lemma 2.2 we can write, 
E((/ M ,/if) 2 ) 



= E ( f P t N f™ ((xi, . . . , *„)) . . . , x 2 „)) ^((fai), • • • , ^(dxan)) 

+ °( — ) 

By Theorem 3.4, 

-> ^/if (0&i, • • • , a^n)) /if ((»n+i, • • • , x 2n )) 

as iV — > oo and is bounded by ||/jf ||oo||/jf \\ooC 2 t 2n ~ 2 where C has been defined in 3.1. 
Hence by the Bounded Convergence Theorem, 

E((/ M /if) 2 ) - / Ptfg «x u ...,x n ))fg ((x n+1 ,...,x 2n )) 

</(0,oo) 2 ™ 

Ho(dxi) . ..n (dx 2n ) 
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as TV — > oo. Now, by Theorem 3.4, 

P t f™((x 1 ,...,x n ))f™((x n+1 ,...,x 2n )) 

= P t f" (( Xl , . . .,X n )) P t fg ((X n+1 , . . .,X 2n )) . 

Hence, as N — > oo, 

E ((/ M , /^f > 2 ) - f / P*/^ • • • , *„)) ^((fai) • • • Mdx n ) 

\</ (0,oo) n 

( / ^t/,? ((^1' • • • ' M^l) • • • M dx n) ) 

\J(0,oo) n J 

= (f M ,^) 2 - 

Step 3: We prove the Proposition 3.3. We have proved before that as M 
{f M ,th)^{f, h)- Hence, 

(/ M ,/i t ) 2 -</,/i t > 2 
as M — > oo. Also, similarly to the way that we proved that 

Urn Urn sup E (|(/(1 - * M (m)),j%)\) = limliminf E (|(/(1 - * M (m)),tf)\) 

1V1 j\j 1V1 iV 

we prove that 

limlimsupE (J(/(l - V M (m)), /if ) = limliminf E (J(/(l - ^ M (m)),tf)f 
Moreover, 

E(\(f* M (m),tf)(f(l-*M(rn)),tf)\) 2 

< E ((/*M(m),^) 2 ) E ((/(l - * M (™))^f > 2 ) . 

Hence, we deduce that 

limlimsupE {\(f* M (m), /if >(/(l - * M (m)), /if )|) 2 

= lim lim inf E (| (/$ M ( m ), /if ) (/(l - ^ M (m)), /if ) |) 2 = 0. 

Hence, we obtain 

limsupE ((/,/if > 2 ) = hminf E ((/,/if ) 2 ) = (/,/i,) 2 . 

AT JV 

Thus, 



hmE((/,/if) 2 ) = (/,/i^ 



as required. 



38 



3.3 Conclusion 

For all r G T and for all / G Cj, (A-(0, t)) we have proved that 

E «/,#»-></, ft) 

and 

E((/,/if) 2 )-(/,^) 2 

as A" — > oo. So we deduce by the remark we did at the beginning of the paper, that for 
all/eC & (A T (0,*)), 

(f,P?)^(f,iH) 
as A" — > oo in probability. Now let us prove our main result. 

4 Proof of Theorem 1.1 

Our aim in this section is to prove our main result stated in Theorem 1.1, that is 

~N 

as A" — > oo weakly in probability. In Section 3, we have proved that for r G T and for all 
/eC b (A T ((M)), 

(/,£?)-></,&) 

as A" — > oo in probability. To prove Theorem 1.1, we need to prove that for all / G 
C b (A(0,t)), 

(f,P?)^(f,iH) 

as A" — > oo in probability. 

4.1 A tightness argument 

We are reviewing here a particular case of the work of [1], [2]. The usual Marcus- 
Lushnikov process gives at each time the distribution in masses of the particles present 
in the system but does not retain any notion of configuration of these particles. We are 
going to define a process (X t ) t > (the Marcus Lushnikov on trees) on T(0, oo), the space 
of trees on (0, oo). The space T(0, oo) is given by 

T(0,oo) = [J T r (0, oo 

where Ti(0, oo) = (0, oo) and for r = {t±, t 2 } G T, 

T T (0, oo) = {y = {yi,y 2 } ■ V\ e T T1 (0,oo),y 2 e T r2 (0,oo)}. 

On T(0, oo) define the mass function w! : T(0, oo) — > (0, oo). For y G Ti(0, oo) = (0, oo), 
we set 

m\y) = y. 
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Recursively for t — {t±, r 2 } G T, for y = {7/1,1/2} G T T (0, 00), we set 

m'(y) = m'(yi) + m'(y 2 )- 

Define also the counting function n' : T(0, 00) — > N as follows. For y G Ti(0, 00) = (0, 00) 
we set 

n'(y) = 1. 

Recursively for r = {t±, t 2 } G T, for y = {1/1,1/2} G T r (0, 00), we set 

n'(y) = n'(yi) +n'(y 2 ). 

Let us define (X t ) t > . Let [N] = {1, . . . , N} and let y 1: . . . , y N > be the masses 
associated to each particle in [N]. Set 

AT 



i=l 



For each i < j G [iV] take an independent random variable Ty such that is exponential 
with parameter K (m'(yi),m'(yj)), and define 



T = mmTij. 

i<j 



Set X t = X for i < T and 

X T = X t - (5 Vi + 5 Vj - 5 {yuyj} ) 



if T = Tij, then begin the construction afresh from X T - Let (X t ) t > be Marcus- 
Lushnikov on T(0, 00) with kernel ^ starting from 

N 

i=l 

Set 

V? = N-'X?. 

Let 

3/ t :A(0,t)->T(0,oo) 

be the map on forgetting times. Hence jxf and Jif are related through the following 
equality 

v? = P?°vr 1 - (22) 

Define 

<p : T(0, 00) — > (0, 00) by (p = ip o m. 
On Ti(0, 00) = (0, 00), <p = <p and so 

(<p 2 ,fi ) = (y? 2 ,^o) < 00 
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by (4). Then [1] tells us that we can find n > such that for all N, 
P (/if ({y G T(0, oo) : n\y) > n }) > C(e)}) < e 
where C(e) > 0. Now by (22) 

/7f ({y G T(0, oo) : n'(y) > n }) = ({£ G A(0,t) : n' (y t (0) > ™o}) ■ 

Hence, 

P (tt e ^(0, : n (y t (0) > n }) > < 6. (23) 

4.2 Proof of Theorem 1.1 

Take f e C h (A(0,t)). For n G N consider 

A n = |J A T (0,t). 

reT 
with ti(t) < n 

Since 

A(0,f) = U^(°'*) 

reT 

it is clear that 

A(0,t) = \jA n . 

n>l 

Moreover, (A n ) n>1 forms an increasing sequence. Set f n = /1a„- Then, f n G C& (A n ) and 
fn — > f as n — > oo. We can write 



/n — ^ /nlA r (0,t)- 

Thus, 



reT 
with n(r) < n 



with n(r) < n 

Now, since f n ^A T (o,t) G C& (A,-(0, £)), by subsection 3.3, 

(fn^A T (0,t), ^ ) — >■ (fn^A T (0,t), fit) 

as A" — > oo, in probability. Thus, since the sum over {reT with n(r) < n} is finite, we 
obtain 

</»,£?>-</»,&> (24) 

as AT — > oo, in probability. Now let us prove our result. Also, for all n, f n — > / as n — > oo 
and |/ n | < |/|. Hence, by the dominated convergence theorem, 

fit(fn) - fit(f) (25) 

as n — > oo. 
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Now let us fix e > 0. Consider, 

I " (f,fit)\ <\(f,tf) - (fn,tf)\ + Kfmtf) - (fn,fh)\ 

+ \{f n ,ih)-{f,ih)\ 

By subsection 4.1, we know that we can find n > 0, such that for all N, 

P (Vf (tt e ^(0, : n' > n }) > < |. (26) 

Now, f-f n = flAo. Thus, 

I (/, ) - </„, ) I = I (/ - fn,fi? ) I < II / II oo | ( U , /if ) | 

For all n > no, C A c nf) . So using the relation (25) we obtain that 

P(K/,£T> - </»,/^>l > C) < P(||/||oo|(U„/if)| > 6) 

e 

for all n > n . Hence, for n > n , for all iV, 

P(l(/^f)"(/n,/if)|>e)<| 
For the third term, by 4.1, we can choose n\ so that, for all n > n 1: 

\{fn,fit) - {f,fit)\ < I 

For the second term : Take n > max(n ,ni). By the relation (24), we can find N so 
that for all N > N , 

V{\(fn,tf)-(fn,iit)\>e)< e - 

Hence, for N > N Q , 

V(\{f,p?)-{f,iH)\>e)<e 

that is 

</,£?>-></,&> 

as N — > oo in probability. Hence, 

~N 

fit fit 

as N — > oo weakly in probability as required. 

5 Appendix 

5.1 A topology on A(0, t) 

We equip A(0,t) with a topology. Define 

4(0,0-^(0,1) 
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as follows. For £ G Ai(0,t) we set £ = 6 Recursively, for r = {ti,t 2 } G T, £ = 
(s, {6,6}) e ^r(0,t) with 6 G A n (0, s) and 6 G A T2 (0, s) we define £ = (f, {6,6}) € 
A r (0, 1) with 6 G A T1 (0, 1) and 6 e A- 2 (0, 1). 

We are going to construct a topology on A(0, 1) and then rescaling we will obtain a 
topology on A(0,t). The set T is countable so we can give it a strict total order <. For 
r = 1 G T, Ax(0, 1) = (0, oo). We equip A^O, 1) with the usual topology on (0, oo). For 
t = {ti, t 2 } G T with T\ < r 2 , we can identify A T (0, 1) to be 

A T (0,l) = (0,l)x A T1 (0,l)x A r2 (0,l) 

Now, for r = {1}, 

A { i } (0,l) = (0,1) x Ai(0,l) x A 1 (0,1) 

We equip (0, 1) with the usual topology on R. We have already given Ai(0, 1) a topology. 
Hence, we equip A{i}(0, 1) with the product topology (Tychonoff topology). By induction, 
assume we have topologies on A n (0,l) and A T2 (0, 1) for ti,t 2 G T. Recursively for 
t = {ti, t 2 } G T with Ti < r 2 , we equip 

A r (0,l) = (0,l)x A Tl (0,l)x A T2 (0,1) 

with the product topology. Then , 

A(0, 1) = (J A T (0, 1) 

is naturally equipped with a topology. Finally, rescaling, we obtain a topology on A(0, t) 
for t > 0. Similarly, we can equip A(0, t) with a topology. 

5.2 A measure on t4|(0,£) 

Fix < i < T. The set T*[n] is countable so we can give it a total order < . Take 
i = {11,12} G T*[n] with i\ < i 2 . Let y,?/ 1 and y 2 by the vector of masses respectively 
associated to %,%i and i 2 . Our aim is to define a measure on j4|(0, i). We are first going 
to construct a measure on Af(0,l) and then by a similar rescalling to the one done in 
subsection 5.1 we will obtain a measure on Af(0,t). Equip Ai(0, 1) with its Borel-u 
algebra £>j. We can We can identify A^(0, 1) to be 

Af(0, 1) = (0,l)x <(0, l)x 4(0,1) 
as ii < 2 2 . For i — 1, for £ G (0, 1), set 

VM) = Sy. 

Recursively for i = {11,12} G T*[n] with ^ < i 2 , with associated vector of masses y,y x 
and y 2 , for | = (s, 6, 6) G A^O, 1), define 

^(deV^(dfi)^(d6)ds. 
This defines a measure on the product space 

B {(0,1)) ®B h ®B i2 

1 2 

where B ix and £> i2 are the respective Borel-a algebra on A\ x (0, 1) and A v i2 (0, 1). Then 
rescalling as in 5.1 we obtain a measure on A?(0,t). 
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5.3 A limit measure on A(0,t) 

Let E = (0,oo). Fix < t < T. Fix n > 0. We define a limit measure on A(0,t) as 
follows. For £ G A k (0,t) with fceN, set 

ji' t (d£) = exp ^- jf J X (y, y'VrW) dr ) Ho{dy) 

where y = m(£) and (/i r )r<T is the strong deterministic solution to the generalized Smolu- 
chowski equation (3). For % = ^2} £ T*N with type r(i) — r — {ti,t 2 } G T , for 
£ — ( s > {6? 6}) e -Ai(0, t), with s < t, define recursively 

VM) = e(r)K(m(^),m(^))il'MiWC2)ew (- J* ^ y V(<V)dr) ds (27) 

where y = m(£) and e(r) = 1 if T\ 7^ r 2 and e(r) = | if Ti = r 2 . 
For £ G A fc (0, i) with fc G N, set 

*€ = 1 

For i = {ii,i 2 } G T*N, £ = (s, {6,6}) e ^(0,t), define recursively 

X e = A-(m(ei),m(6))^ € i^6 
Theorem 5.1. Let i = {h, i 2 } G T^N with type r = r(i) G T. 

1. TTien /or £ = (s, G ^(0,*), 

= 2~ q{T) Kz exp ( - / / K(y r , y')fi U{r) (dy')dr J i/f (d£)n (dyi) . . . fM)(dy n ) 
\ Ja(0 Je J 

where y — (y±, . . . , y n ) is the vector of masses associated to £ and uf is the measure 
described in subsection 5.2. 

2. Take f G Cb(A T (0,t). Let fi = fog i where g,i : Ai(0,t) — > A T (0,t) is the map on 
forgetting labels. Then, 

[ minm = [ m) [ mikm 

JA T (0,t) JAi(p,t) JA T (0,t) 

Proof of Theorem 5.1: 

(a) Let us do it by induction. It is clearly true for % G T*N. Fix n > 0. Suppose it 
is true for all % G TJlN with k < n — 1. Is is true for k — nl Take i = {ii,i%} G T*N. 
Then, writing ra(ii) = and n(i 2 ) = n — k we have 

n(ii) < n — 1 

and 

n(? 2 ) < n — 1 
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Take £ = (s, {£1, £ 2 }) £ A(0, t). Without loss of generality, assume that %\ is formed from 
the particle 1, . . . , k with associated masses y l = (yi, . . . , y^) and that i% is formed from 
the particle k + 1, . . . , n with associated masses y 2 = (yk+i, ■ ■ ■ , y n )- Using the induction 
hypothesis, we have, 

fi'Mi) =2- g(Tl) % exp (- f [ K(y r ,y>) mr) (dy')dr 

\ J A(£i) JE 



»i 

and 



A(6) 

^ KiW^i) ■ ..(io(dy k ) 



ji' s {d&) =2-^ 2 ^ 2 exp f- / / K(y r ,y')nn(r)(dy')di 

V ./A(£ 2 ) •/£ 



'A(6) 

2 

i£ (d&)no(dy k +i) ■ ..Ho(dy n ) 



Using the relation (27), we have 

#(d0 =c(r)2-^)2-^)^Mei),m(6))/Q 1 ^ & i/if (dei)i^ • • -Mofc) 

exp ( - / / K(y r ,y')ii u{r) (dy')dr) exp ( - / / K(y r ,y')fj, n{r) (dy')dr 
V ^A(fi) /V JA(6) 



exp ( - / / K(y,y')fj, r (dy')dr 

Js JE 



Now, K 5 = X(m(6),m(£ 2 ))X a ^ 2 , z/fK) = ^(dfc)^ A(f) = A^UA^U 

((s,t] x f) and 2~^ T ) = 2" 9 ( ri )2-^ r2 )2- e ( r ). Hence, 

= 2"«W^ exp f- / / X(y r , j/'Km^')^ ■ ■ ■ Ho{dy n ) 

\ Ja(o Je J 

as required. 

□ 



45 



5.4 Some simulations 



Here are some simulations under Visual Basic of the Marcus-Lushnikov process on trees. 
The graphics below represent trees that have been simulated following the Marcus- 
Lushnikov process on trees with different kernel K and an initial number of particles 
N. In these simulations, all the initial particles have mass 1. These pictures show for 
each kernel the sort of trees limit we can expect to find in the limit measure. 




Figure 1: N = 128, K(x,y) = 1 





Figure 2: N = 128, K(x,y) = xy 
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Figure 3: N = 128, K(x, y) = l/(x + y + l) 
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